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Abstract 

We study the semi-classical trace formula at a critical energy level for a /i-pseudo- 
differential operator whose principal symbol has a unique non-degenerate critical 
point for that energy. This leads to the study of Hamiltonian systems near equi- 
librium and near the non-zero periods of the linearized flow. The contributions of 
these periods to the trace formula are expressed in terms of degenerate oscillatory 
integrals. The new results obtained are formulated in terms of the geometry of the 
energy surface and the classical dynamics on this surface. 
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1 Introduction 



Let be Ph a /i-pseudodifferential, or more generally /i-admissible (see [15]), self- 
adjoint operator on M". Tlie semi-classical trace formula studies the asymp- 
totic behavior, as h tends to 0, of the sums 

^{E,h,^)= E ^i^^^), (1) 

\Xj{h)-E\<e 

where the Xj{h) are the eigenvalues of Ph- Here we suppose that the spectrum is 
discrete in [E — e, E + 6], some sufficient conditions for this will be given below. 
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Let be p the principal symbol of Ph and the Hamilton flow of p. The semi- 
classical trace formula establishes a link between the asymptotic behavior of 
(1), as /i — >■ 0, and the closed trajectories of $t of energy E. An energy E is said 
to be regular when Vp(x,^) 7^ on E^, where = {{x,^) / p{x,^) = E} 
is the surface of energy level E, and critical if it is not regular. The case of a 
regular energy has been intensively studied and explicit expressions in term 
of are known for the leading term of (1), under suitable conditions on the 
flow and when the Fourier transform of is supported near a period, see, 
e.g., Gutzwiller [7], Balian and Bloch [1] for the physical literature, and from 
a mathematical point of view Brummelhuis and Uribe [3] , Petkov and Popov 
[14], Charbonnel and Popov [4], Paul and Uribe [13]. 

Here we are interested in the case of a critical energy Ec of p. Brummelhuis, 
Paul and Uribe in [2] have studied the semi-classical trace formula at a critical 
energy for quite general operators but limited to "small times", that is for 
supp((^) contained in such a small neighborhood of the origin that the only 
period of the linearized flow in supp((^) is 0. Khuat-Duy, in [10], [11], has 
obtained the contributions of the non-zero periods of the linearized flow for 
arbitrary (f with compactly supported (p, for Schrodinger operators —A+V{x) 
with V{x) a non-degenerate potential. In this case the main contribution of 
such a period was obtained as a regularization of the Duistermaat-Guillemin 
density pt{x,^) = \dct{d^t{x,^) — ld|^2. Generalizing Khuat-Duy's result to 
more general operators was an open problem and is the purpose of this article. 
For a critical energy level Ec of an arbitrary /i-admissible operator, and (a;o, Co) 
a critical point of p, the closed trajectories of the linearized flow 



are not necessarily generated by a positive quadratic form, contrary to the 
case of a Schrodinger operator, and will give rise to new contributions to the 
semi-classical trace formula, other than those obtained in [2], [10] and [11]. 
More precisely, viewing ^^.^^^^(xo, {0) as the Hamiltonian flow of the Hessian 
d^p{xo,^o), here interpreted as an intrinsic quadratic form, these new contri- 
butions to the trace formula arise from the non-trivial closed trajectories of 
4,e*t(^o,Co) of zero energy 



The reader can easily verify that the set of such u is empty in the case of 
a Schrodinger operator : this explains why the new contributions obtained 

here does not appear for these operators. We will show that the new contri- 
butions are supported in {Ec,Xo,C,o,T,u) with (T,u) satisfying (3) and can 
be expressed in term of d'^p{xQ, ^0) and higher order derivatives of the flow in 



<Pt{u) = 4,«$tK,eo)« = u, ue T(,„,5„)(T*M"), t 0, 



(2) 




Mu) =u,T^ 0, 
■p{xo,^o){u) = 0, uj^Q. 



(3) 
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2 General hypotheses and main results 

Let Ph be in the class of /i-admissible operators on with a real symbol. We 
refer to [15] for the principal notions of semi-classical analysis which we will 

use. We note p the principal symbol of Ph, and the sub-principal symbol. 
For Ec a critical energy of p we will study the asymptotic behavior of the 
spectral function ^{Ec, h), defined by 

^{E,,h)^ E (4) 

Xj{h)e[Ec-e,Ec+e] 

under the following hypotheses which are classical in this context : 

(Hi) 3^0 > such that p^^{[Ec — Eo, E^. + £o]) is compact. 
{H2) zo = {xq,^q) is the unique critical point of p on the energy surface ^Ec- 
{H3) zq is a nan degenerate critical point of p and its Hessian (i^p(^o) is 
diagonal in some suitable set of local symplectic coordinates near zq 

p(x,0 = ^c + ^Eti^,((^,-^oj)'+<7,(0-eoj)')+C?(||(x-Xo,e-eo)||'), (5) 
i=i 

with (Tj — ±1 and wj e R\{0}. 

{H4) (f is in the Schwartz space S{M.) with Fourier transform (p e C^(R). 

By a classical result, see e.g. [15], the hypothesis (i/i) insures that the spec- 
trum of Ph is discrete in — \Ec — e,Ec ^ e\ for £ < £0 and h small 
enough : this will be assumed in the following. We note Exp(tiJ/), with 
= d^f.d^ - d^f.d^, the Hamihon flow of a function / e C°^(T*R"). For 
$t = Exp{tHp) taking the derivative with respect to the initial conditions 
gives a symplectomorphism d^^^^tix,^) : T(a,^^)(T*M") T$,(a,_g)(T*M"), and 
for Zo — {xo, ^0) a critical point of p we have the fundamental automorphism 

4,€*t(-2o) : 7;„(r*M") ^ T,,{T*R^). (6) 

Near the critical point of p we can write 

N 

P{X, = + E PJ i^,0 + 0{\\{x-Xo,^-CoW^'), 

where the functions pj are homogeneous of degree j in {x — Xq,^ — ^q). In 
particular, p2 is the Hessian in zq and can be interpreted as an invariantly 
defined quadratic form on T^^{T*W). 

Definition 1 For all T e R let be = Ker(4,^$T(;2o) - Id) C T,^{T*K^) 
and — T'2;g(T*R"')/5'r- To this linear subspace we associate its dimension 
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It — — dim(5'T); and also the following three objects 



CQ, = {(x,Oe^T/QT(x,0 = 0}. 



2?T =P2|^r> 



(7) 
(8) 
(9) 



We say that T is a period of dx,^^t{zo) if dr 7^ {0} and that T is a total 
period ofd^,^^tizo) when dr = T,,{T*W). 

The next condition is inspired by proposition 2.1 of D. Khuat Duy [10] con- 
cerning Schrodinger operators and will be useful to separate the contributions 
of fixed points from those of the non-trivial periodic trajectories. 

(i^s) For all period T of d^^^^ti^o) there exists neighborhoods Vt of T and 
Ut of Zq such that 7^ z for all T & Vt and all z e Ut\{zo} H E^^. 

Note that (H^) only concerns the dynamics on the energy surface J^e^ and 
Khuat Duy loc. cit. has shown that (ifs) is always satisfied by a Schrodinger- 
type Hamiltonian -|- V{x). We give in the last section two examples of non- 
Schrodinger operators satisfying (H^), to show that this class is non-empty. 
We note a the usual symplectic form on T2„(T*]R") and by the derivative 
of order /, with respect to initial conditions, evaluated in {t,z). If m is in a 
vector space V we use the notation = {u, ...,u) e VK 

Definition 2 If k > 1 is the first integer such that d^^T 7^ 0, we put 



sphere, where we are working in local coordinates such that (5) holds. Finally, 
let dLq^ be the Liouville-measure on this surface, i.e. dLq^ A dQT = dO on 
Cqj. n §2dT-i^ where dO is the surface measure ofS'^^'^~^. 

The derivatives 6?^^$; are computed in section 4.4 and the relation (i^g$r 7^ 
can be stated in terms of resonances, i.e. the arithmetical properties of the 
eigenvalues of p2 or, more precisely, of Qt- 

If © is a cut-off function near the energy Ec we have 



(10) 



and 




Ph- 
h 



)e(p,)) + 0(O, 



(12) 
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as will show proposition 6 below. Hence, modulo 0{h°°), we can write 

7(£;c, h) = Tr(C¥'(^^^)e(PO) + ^((1 - KM^^^MPk)) 

= -fi{zo, Ec, h) + 72(2:0, Ec, h), 

where is a function with compact support near zq such that {H^) is valid 
on supp(V') and the notation t/j'^ stands for the classical Weyl /i-quantization. 
Under the additional hypothesis of having a clean flow, the asymptotics of 
j2{zo, Ec, h) is given by the regular trace formula. 

Without any loss of generality we can suppose that supp(<^) is small enough 
near some non-zero period T of d^t{zQ) such that T is the only period of 

d^t{zo) on supp((^), as follows by an easy partition of unity argument, the set 
of periods of the linearized flow being discrete. For the remaining contribution 
71(2:0, Eci h) we then obtain 

Theorem 3 Under {Hi) to (H^) and the assumption that Qt is positive or 
negative definite, we have 

71(^0, Ec, h) = C{T)Kt{v) + 0{h^, ash^O, 

where 

n(rr\ 1 exp(2f sgn(gT)) dr - 1 . , . , . 

C{T) = -- ] 1 exp(^7r— — sign(gr))r(cir), 

^ |det(gT)h 

and 

[2Tiy+^T \ |det(Id-(i$t(2:o)|2 / 



In the case of a non definite Qt we have : 

Theorem 4 If n > 2, under (Hi) to (Hr,) and if Rk{z) ^ for all z e 
Cqj,\{0}, we have 

71(^0, Ec, h) = h— ''^{Kt0{T) exp(^V(zo)) + 0{h-^)), for h ^ 0, 

where 

drp — 1 ~ f ~ 2drj.-2 

KT = fJ.k{T)exp{i7r^j^sign{Rk)) J |i?fe(^)| —dLQ^{9), 

and 



pik{T) = --r 



1^^2rfr-2 exp(ifsgn(gT)) 



k k ' \det{qT)\H2TrYT+i 
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More generally, if {z e Cq^\{0} : Rk{z) = 0} ^ 0, but if VQT{z),VRk{z) 
are linearly independent on this set, then the same result holds with 



Remark 5 Theorem 3 contains the class of Schrodingcr operators + V{x), 
with V non degenerate, since in this case p2|Ker((i$t(zo) — Id) is positive 
definite. Moreover, the condition n > 2 in Theorem 4 is natural, since, under 
om^ assumptions, in dimension n = 1 we have P2{,xi,^i) = ^Wi{xf + $^f), or 
^Wi{xl — ^f). In the first case Qt is definite and in the second case is the 
only period of the linearized fiow. 



3 Oscillatory representation of ^{Ec, h) 



We introduce a cut-off function O G C^{\Ec — e,Ec + e[), such that = 1 
near Ec and < G < 1 on M. We then have the decomposition 



^i(E,,h)+y'{E,,h). (13) 



Proposition 6 YiEc, h) = 0{h°°), when h ^ 0. 

Proof. With y? G S(R) for all A; G N there exist Ck such that \x'^(p{x)\ < Ck on 
M. If N{h) is the number of eigenvalues inside [Ec — e, Ec + s]r\ supp(l — 0), 
then, by Theorem 3.13 of [15], we have the estimate 

\j'{E,,h)\ < N{h)C,\ ^'^''^~ ^' \-^ N{h)^0{h-^). 

On the support of (1 — 6) we have \Xj{h) — Ed > Eq > 0, this gives 

W{Ec,h)\ < N{h)CNe^^h^ < cjv/i''"". 
Since it is true for all /c G N the result follows. ■ 
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As a consequence, the asymptotics of 'y{Ec,h), modulo 0{h°°), is given by 
^"{Ec,h). Now, by Fourier transform, and inversion, we have that 

R 

The trace of the left-hand side is then exactly j"{Ec, h), and we can write 

Y{Ec, h) = ^Tr j Q{Pn)ExpC-lPh)e-''^mdt. (14) 

The formula (14) uses the localized unitary group UQ^h{t) = 0(P/i)Exp(^P/j). 
A classical result, see e.g. [5], about this object is 

Proposition 7 Let A he the Lagrangian manifold associated to the flow of p 
A = {(t, r, X, e, I/, 11) e 7^*^'"""' / r = -p(x, 0, {x, i) - ^t{y, r?)}, (15) 

then U@^h{t) is a h-Fourier integral operator (or h-FIO) associated to A. More 
precisely, there exists for each N a FIO UQ^fl{t) with integral kernel in the 

Hormander class /(M x x ]R",A) and operators Rl^\t) on L^(R"), with 
uniformly bounded norms for < h < 1 and t in a compact subset of M., such 
that 

UeAt)-U^e7(t) + h''Ri''\t). (16) 

We recall the compactly supported cut-off function ip = iIj{x,C,), on whose 
support (i/5) holds. If ipi is such that ipiip = ip, with supp('^i) small enough, 
then by cyclicity of the trace, 

72(^0, ^c, h) = ^Tr J ^(i)Ce(P,.)Exp(|(P, - E,))i;r,^dt + 0(h^), 

and i/j^e{Ph)Exp{'^{Ph - £^c))^]^ft e /(M x R" x M", A). After perhaps a local 
change of variable in y, the operator ■0^©(P^)Exp(^P^)'05^;j can be approxi- 
mated , modulo an error 0(h^), by an h-FIO with kernel 

I em'^^^^^-^y'iHr,{t,x,y,^,h)d^, 

sec 4.1 below. Integrating this kernel on the diagonal gives, modulo terms of 
order 0{h^) 

R T*R" 

aN{t, X, ^, h) = (f{t)bN{t, X, X, ^, h). (18) 



7 



For a detailed construction we refer to [15] or [5]. 

Remark 8 Because of the presence of Q{Ph) and ip, the amphtudes gn are 
of compact support. Since we are interested in the main contribution to the 
trace formula we note a the amplitude of (17), i.e. a depends on {h, N). 



4 Study of the phase function and of the classical dynamics 



First, we study the nature of the critical points of (17). After we establish 
some results on the classical dynamics related to p{x,^) and we compute 
the Taylor expansion of the phase function. Resonance-type conditions will 
naturally occur in the study of this question. 



4-1 Singularity of the phase 



By Theorem 5.3 of [9], we can, after perhaps a local change of variables in y, 
suppose that the flow near {xq, ^o) and for t e supp(<^) sufficiently small, 
has a generating function S{t, x, r)), i.e. 



{x,c)^My,v)^ < 



y = d^S{t,x,ri), 
^ = d^S{t,x,r]), 



(19) 



which therefore, by a classical result, satisfies the Hamilton- Jacobi equation 
dtS{t, X, rj) +p{x, dxS{t, x, rj) = 0. Hence, near (T, xq, ^q), the Lagrangian man- 
ifold A of the flow is parameterized by the phase function S{t,x,r]) — {y,r]). 
This choice for the phase is only valid near (xq, Co) when .^o 7^ 0, but if = 
we can change the operator Ph by e^^^'^^^P/^e^'^^^'^i^ with ,^1 ^ 0. This does 
not affect the spectrum and the new operator obtained has symbol p{x, ^ — d) 
and critical point (a;o,Ci)- With (19) a critical point of (17) satisfies 



d,s{t,x,O^C, 

d^S{t, x,^) = x, -v^ 
dtS{t,x,0 = -Ec. 



$,(x,0 = (x,0, 

p{x, = ^c- 



The following lemma is classical and can for example be found in [11]. Recall 
that we will often denote points (x,^) of phase space by a single letter z. 

Lemma 9 Let us define ^'(t, a;,^) = S{t,x,C,) — (x,^) + tE^, then if Zq is 
critical point of "if, we have the equivalence d1'if{zo)Sz = 4^ dz^t{zo)Sz — 5z, 
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y5z G T2(,(T*]R"), i.e. the degenerate directions of the phase correspond to fixed 
points of the linearized flow at Zq. 

If we use Lemma 9, we obtain for our phase function 

Corollary 10 A critical point (T, a;o,^o) ^) is degenerate with re- 

spect to (x, ^) if and only if T is a period of the linearized flow (ia;,c^t(2^o, Co)- 

The next result is also well known from classical mechanics 

Lemma 11 If dxp{xQ,^Q) = d^p{xQ,C^Q) = then (i^^g$t(xo, Co) is the Hamil- 
tonian flow of the quadratic form |<i^p(a;o, Co) on Txq,^q{T*M.'^) . 

Let T 7^ be a period of d$t(zo). Corollary 10 shows that we must introduce 
X, = S{t, X, C) - {x, 0-{t- T)g{t, X, C) + {S{T, x, C) - {x, C)). (20) 
This function g is C°° and satisfies 

(99 

g{T,x,0 = -g^iT,x,0 = -p{x,d,S{T,x,0). 
To simplify the notations we write g{t, z) = g{t, x, C) and 

R{z) = R{x, C) = S{T, X, - {x, . (21) 



Lemma 12 In a neighborhood of zq, and near T, the only critical point, on 
the energy surface T^Ec of the functions S{T, x, C) — {x, C) and g{t, x, C) is zq. 

Proof. First, dx,^{S(T, x,^) — {x, ^)) = is equivalent to ^t{x, C) = {x, C)- But 
near zo and with (H^) this can only be satisfied for (x,^) = zq. 
Next, we consider g{T,x,^) = —p{x,dxS{T,x,^)). here dx,^g{T,x,^) = is 
equivalent to 

a,p(x, d,S{T, X, 0) + d^p{x, d,S{T, X, mdl^T, X, C)) = 
d^p(x,d,S(T,x,mdl^SiT,x,0) = 0. 



Clearly zo is critical since dxS{T, xq, Co) — Co- If is a critical point of g{T, x, C) 
we have 



i„e^5(r,x,c) 



1^0 dl^S{T,x,i) 



dXx,d,S{T,x,C)) 
d^p{x,d,S{T,x,^)) 



0. 



(22) 
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Where I„ is the identity matrix of order n. But ^tid^S{t, x, C),0 — (^^ dxS{t, x, ^)) 
and this leads to 



d,,^Md^S{t,x,C),C) 



^dl,S{t,x,0dl^S{t,x,0\ (in ^ 



V In J \dl,Sit,x,Odl^Sit,x,OJ 

Since d^t is an isomorphism Eq.(22) imposes 

d,p{x, d,S{T, X, 0) = d^p{x, d,S{T, X, 0) = 0. 

In a suitable neighborhood of Zq this imphes that (x,^) = zq, since zq is non- 
degenerate. With {t — T)g{t,z) = S{t,z) — S{T,z), for t T, the critical 
points of g{t, z) are those of S{t, z) — S{T, z), and equation (19) shows that 

Md^S{t, X, 0,0 = i^, d.Sit, X, 0) = {x, d^SiT, X, 
= ^T{d^S{T, X, 0,0 = Md^S{T, X, 0,0- 

We then have by the group law 

^^tiMd^SiT, X, 0,0) = ^T-t{d^S{t, X, 0,0 = (^^-^l^, ^, 0,0)- 

For t ^ T and |t — T| smaU the point {d^S{t, a;, 0, 0) would be periodic, with 
period {t — T). For (x, ^ ^Ec) near 2:0, {H^) implies that (a;, = zq. ■ 

The Hessian matrix with respect to 2; = (x, of g in (T, 2:0) satisfies 



(Hess,(y)(T, z^){5x, 50, (5^, ^0) = (Hess(p)(;so)i?(5x, 50, ^(^x, ^0) , 
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In 



(23) 



dl.S{T,z,)dl^S{T,z,)) 
with S non singular, as seen before. This proves the relation 

Hess,(^)(r,^o)= *5Hess(p)(^o)5- (24) 
In the case of a total period T of d^t^^o) we obtain 

Proposition 13 Ifd^tizo) is totally periodic, with periodT , then the function 
g{t,x,^) satisfies }iesSz{g){T, zq) — Hess(p)(2;o). 
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4-2 The linearized flow in zq 



Up to a permutation of coordinates we can assume that 

The flow of p2, viewed as an element of End(ro(r*R")) ~ End(R2"), is 

^ a{t) c{t) ^ 

a; 



(25) 



Y.MtH,,){x,i)=A{t) 1^1, 



h{t) d(t) 
-c(i) a{t) 
y h{t) I 



where (x, = (2;', x", C', i"), x' , i' e M^ x" , C" e 



pn— A; 



, and 



a(t) = diag(cos(wit)), b(t) = diag(ch(wit)), 
c{t) — diag(sin(wji)), d{t) — diag(sh(wii)), 



and " diag" means diagonal matrix. In the following we work on the subspace 
{x" = ^" = 0} obtained by projection on the periodic variables. Let be / a 
subset of {1, n} with / elements, / > 1. The existence of a non trivial closed 
trajectory of dimension / imposes that there exists a c e M* such that 



Vi e /, 3ni e (Z*)' : Wi — crii. 



(26) 



Remark 14 Let M{w) = {k G / {k,w) = 0} be the Z-module of reso- 
nances of the vector w = {wi, ...,Wn)- The relations (26), for / > 1 lead to reso- 
nances, since riiWi — rijWj = 0, but a resonant system can have no periodic tra- 
jectories of dimension greater than 1, as is shown by {^/2 + ^/3, V2 — ^/3, ^/2). 

Since we are interested by the periods of d^ti^o) we define (assuming (25)) 

1 ^ 

and also 

in- 

^ \w\ 
i=fci+i ^ 
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so that Qper{x,^) = Q~^{x, C)~Q i^, 0' Lct Pi and P2 be the hnear subspaces 
obtained by projecting orthogonally on the effective variables of Q~^, Q~. 

Proposition 15 If we have C Pi or 5t C P2, then }iess{g){T, zo)\:sj. is 
respectively positive or negative definite. 

Proof. Choosing coordinates as in (25), we see that for v & 



^^^^^dl,S{T,zo)dlS{T,zo) 
L, 



\dl.S{T,zo)dl^S{T,zo)) 



V. 



With dim(g^r) = 20?^ we can choose our coordinates such that 



a(T) , ^ 
b{T) ' 







^0 * 
V b{T) 



Elementary considerations show that 



dl^SiT, zo) 



* 



* 



(27) 



where * designs matrix blocs which are irrelevant for the present discussion. 
Hence, by restriction to '■ {^^ss(g)(T, zo))\:g^ — (Hess(p)(2;o))|5T- * 



4.3 Taylor series of the flow near Zq 



We start by the general case of an autonomous system near an equilibrium. Let 
be $t the flow of a C°° vector field X on R" with coordinates z — {zi,...,Zn) 
and let zq a fixed point of We denote by ^(^0) the matrix of the lineariza- 
tion of X in Zq 

dX^ 

^(^0) = (^)(ao(^o), (28) 

and we recall that d^t{zo) = Exp{tA{zo)). Here, and in the following, the 
derivatives d will be taken with respect to z, we denote by rf^/ the k-th 
derivative of / regarded as a multi-linear form on the fc-fold product MJ^ x 
... X IR'*, and d'^f{zo) or d'^^f this derivative evaluated in zq. As an example. 
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we compute the second derivative of the flow in zq : 



'dt ^^dzkdzi dzj dzi 

Hence at the point Zq we obtain 

Let us write Hcss(X)(2;o) for the vector valued Hessian of X evaluated in zq. 
Interpreting (29) as an inhomogeneous system of equations for ^z^,zj^t{z^ we 
obtain, since $0 — Id and therefore ^'^ J^^^iz^ = 0, that 

t 

j2 



d^^t{zo){z,z) =d^t{zo) I d^_,{zo)}iess{X){zo){d^s{zo){z),d^s{zo){z))ds. 

(30) 







Now, let us assume that Zq is the origin, we generalize as follows : 

Proposition 16 If X{z) = Az + Xk^z) + C>(| |^| 1*^+^), where Xk{z) is homo- 
geneous of degree k > 2 in z, then the flow of X satisfies 

d^MO)^0, yje{2,...,k-l}, 

t 

d'^t{0){z,...,z)^d^t{0) J d^^s{0)d'X{0){d^s{0){z),...,d^s{0)iz))ds. 



Proof. The first result is trivial. At the order k we have for |q;| = A; 

Hence, for 2; = we simply have 

|((^)"*tW)|.=o = c^"^(0)(c?$*(^o),...,ci**(^o))+^((^)"**(^))|.=o. (31) 

Eq. (31) is linear of the form dtUa{t) = faif) + ^^^(t), and by integration, 
with the initial condition ^^=$0(0) = 0, we obtain 



^;)"<l>t)(0) = C(<l>t(0) 



and the result holds by linearity. 
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A more general result is 



Theorem 17 Let be Zq an equilibrium point of X and $i the flow of X. For 
all TO e N*; there exists a polynomial map Pm, of degree at most m, such that 





(32) 

In addition is uniquely determined by the m-jet of X in Zq. 

Proof. For m = 1, d^t{zo) is determined by the operator A{zo), i.e. by the 
1-jet of X. We note a;' G (M")' the image of x under the diagonal mapping, 
with the same convention for any vector. If / and g are smooth we obtain 

m 

d^{fog){xo){x^) = E E c^dif{g{x^mdg{x^){x)rn^...^ {d^g{x,){x^)r^rn)^ 

3=1 aelj 

m 

with Ca e N, Ij = {a e / J2 kaj^ = j}. Since zq is a fixed point, we have 

k=l 

rn 

d^{Xo^t){zo){zn = E E Cad^X{zo){{d^s{zo){z)r^,...,{d^^s{zo){znp-). 

j=l aelj 

For Y — (Yi,..., Ym), we can define 

m 

Pm{Y) = E E Cad^X{zo){Y,"'\....,Y;^-)-dX{zo){Yrr,), (33) 

j=l aelj 

this leads to the differential equation, operator valued 

^(d"'$t(^o))(0 =^(^o)c^™*t(^o)(^™)+Pm(c^$s(^o)(^),...,(i'"-'*s(;2o)(;2™-')). 

With the initial condition d"*$o(^o) = 0, we obtain that the solution is given 
by (32). Moreover, Eq.(33) shows that Pm is completely determined by the 
derivatives of order less or equal than to of X. ■ 



4-4 Application to Hamiltonian systems 

Proposition 16 applied to the fiow of Hp shows that 

t 

d'^^t{zo){z,z) = d^t{zQ) / d^-s{zo)^fiss{Hp){zQ){d^s{zQ){z),d^s{zQ){z))ds. 
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Wc now consider more closely (P^tizo) for t — T, a period of d^tizo)- We 
introduce the following terminology 

Definition 18 ty e (M*)"- is pseudo-resonant to the order I eN if 

I 

3 {ii, ...,ii) e {1, ...,n}, {si^, .. e {-1, 1} such that '^£ijU)i^ = 0. (34) 

Remcirk 19 This notion is weaker than the usual resonance condition since 
for I even there always exists a pseudo-resonance. For example, for Z = 4 we 
have {wi — wi) ± {wj — wj) = 0, although w can be non-resonant at the order 
4. We also observe that all resonances of order 3 are pseudo-resonances. 

In term of the frequencies Wi, we then have 

Theorem 20 // the frequencies w satisfy no pseudo-resonance relation of or- 
der 3 and if T is a total period of d^ti^o), we have (i^$r(^o) = 0. 

Proof. Under the condition (H^), d'^^ti^o) can be expressed as a linear combi- 
nation of integrals of the elementary functions s i— exp{±is{wi + eiWj-\-e2Wk)) 
with Sj = ±1. Hence, to determine (i^$T(^o) we must compute 

T 

j exp(±is(wi SiWj -\- £2Wk))ds, 



but under the assumptions of the theorem all these integrals are 0. ■ 

For all / G let be Mi{W) = {k e Z^" / {k, {w,w)) = 0, |A;| = /} the Z- 
module of resonances of order /. In the presence of resonances we can say the 
following : 

Proposition 21 For k e M^iW) let be k ^ {\ki\,..., \k2n\)- If^k e M^iW) 
d^p 

we have — ^{zq) — and if T is a total period of dz^^t then d'^^Ti^o) — 0. 
dz'' 

The proof is trivial when going back to the proof of Theorem 20. ■ 
Let /* be the puUback by a map /. Then by proposition 16 we have : 

Corollary 22 For a Hamiltonian system with the equilibrium point Zq and 
such that d{^p = 0, Vj e {3, ...,k — 1}, we obtain 

d'^t{zo)^o, yt, yj e{2,...,k-2}, 

t 

d'-'Mz'-') = d^tizo) J d^_s{zo) (d^s{zond1;'H,)) {z''-')ds. 
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And Theorem 20 generalizes trivially to the order k under the conditions of 
Corollary 22. More precisely, if k is odd and if there is no pseudo-resonance of 
order k then, under the assumptions of Corollary 22, we have d^~^^t = 0. 

4-5 Relation between the phase and the flow 

Like in the preceding section we consider a Hamiltonian function p with total 
period T for the linearized flow, satisfying, near 

p{z)^E, + p,{z) + 0{\\z\\'^). (35) 

We recall that z = {x, ^) and z^ = {z, z) e M^"*^. By Taylor, we have 

{k-l)\ 

Under these conditions we can write the generating function at time T as 

5(T, ,T, = {x, + Rk{x, + Rk+i{x, 0, (36) 

where Rk is homogeneous of degree k and Rk+i is the remainder of the Taylor 
expansion. Let J be the matrix of the standard symplectic form on T*M'^. The 
relation between the phase function and the flow of p is given by 

Proposition 23 Under conditions (35) the {k-l)-st derivative of the flow <^t, 
at time t — T , equals : 

d^-i$T(0)((x,O'-') = -{k - l)\JVRu{x,0 (37) 
In addition we have Rj{x, ^) = 0, 3 < j < k, and 

1 ^ 

Rk{x, = ^ / ^((^, 0, d^-siO) {d^.iOYidt'H,)) {x, 0'-')ds. (38) 
■ 

Proof. With equation (19) we obtain 

^^{x+d^Rk+d^Rk+uO = i^+diRk,0+jj^^d'-'^Tm^^^ 
by identification of homogeneous terms we have successively 

{-d^Rk,d,Rk){x,0 = -JVRk{x,0 = jj^^d'-'^Tm{x,C)''-'), 

Rk{x,0 = I {{x,0,-J^Rk{x,0) = l^c7{{x,o,d''-'M0){{x,o'-')), 

where the last result holds by homogeneity. Corollary 22 then imphes (38). ■ 
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5 Normal forms of the phase function 

In this section we derive suitable normal forms for ^{t, x,^) = S{t, x, ^) — (x, ^) 
in our oscillatory integral representation of 72 (-Ec h). We recall the decompo- 
sition 

^{t,x,0 = {t - T)g{t,x,0 + R{x,0, 
cf. formulas (20) and (21). In the micro-local neighborhood of Zq = {xq,C,o) we 
are interested in, the only critical point of R and of g{t,-) is zo for t close to 
T and, moreover, zq is non-degenerate for the latter. 

A further very important simplifying assumption we will make for the moment 
is that, until further notice, T is a total period of d^ti^o). We will show in 
section 6.3 below how to relax this assumption. If T is such a total period, 
then clearly R{z) — 0{\\z\\^). This can be made more precise 

Lemma 24 // near Zq the function p satisfies (H-^) and condition (35) then 
for t near T there exist a non degenerate quadratic form Qt{x,^) such that 
Qt{x,0 =P2ix,0 and 

S{t, X, - {x, 0^{t- T) {Qt{x, + h{t, X, 0) + R{x, 0, 

withR{x,i) = C(||(a;,OII*^) andh{t,x,C) = C(||(a;,^)||'^) uniformly int. 

Proof. On replacing t — T by t we can write \[' = R{z) +tG{t, z) with G{t, z) = 
g{t+T, 2;). By a second order Taylor expansion around zq and proposition 13 we 
have that Glt,z) = Qt{z) + h{t, z), with Qo{z) = P2{z) and h{t,z) = 0{\\z\f). 
Now by proposition 23, R{z) — C(||2;||'^), given that p satisfies (35), and since 
^r^(z) = z + 0(11^11'=-^), we have that S{T,x,C) = {x,^ + 0{\\{x,C)\\''). 
Therefore h{t, z) — 0{\\z\\^) uniformly in t and the lemma follows. ■ 

Reduction of the phcise with respect to Cq^,. 

Let us suppose S'(T, x, ^) contains effectively some terms of order k. We write, 
as before, R{z) = Rk{z) + Rk^i{z), where R^ is the homogeneous component 
of degree k oi R and -R^+i is the remainder of the Taylor series. The following 
lemma is useful for any perturbation by a function of odd degree. 

Lemma 25 Let be Q a non degenerate quadratic form on W^, n > 3, with 
inertia indices greater than 2. For all odd continuous function R, Q and R 
have a common zero on S""-*^. 

Proof. Up to a linear change of coordinates we can assume that Q{x) = \ |^— 
||a;2|p, with Xi G M.^, X2 G M'^, p,q > 2, p + q = n. Now the cone of the 
zeros of Q is invariant under isometrics of the subspaces (a;i,0) and {0,X2). 
By rotating around the origins there exist a continuous curve 71 inside the 
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cone of Q mapping {xi,X2) to {—x\,X2) and a curve 72 mapping {—xi^x^) 
to (— —3:2), inside the cone. If 7 = 71.72 is the union of the two previous 
curves, the function i?(7) gives the result by continuity since R is odd. ■ 

Remcirk 26 A consequence of Lemma 25 is that the set Cq^ n Cr^ fl S"~^ is 
not empty when the function is non-zero and odd and Qt is non-definite. 

We choose polar coordinates z = {x,$^) = r6, 6 G S^"^^(]R). These coordinates 
will perform a "blow-up" of Mx (T*M"\{0}). In general one uses the projective 
space P2„_i(M), but here, since the singularities are carried by the conic set 
of the zeros of Qt, it is convenient to use the sphere §^"~^. For any function 
/, positively homogeneous on M", we note C/ = {x G / f{x) = 0}, the 
conic set of the zeros of /. Finally, g — h means that applications g and h are 
conjugated by a local diffeomorphism. 

Lemma 27 If 9o G S^"~^ with 9q ^ Cq^, there exists a system of local coor- 
dinates X, near {t,r,9) — (0,0, ^0), such that ^ ~ XoXi ^ neighborhood of 
(Xo,Xi) = (0,0). 

Proof. Using the notations of Lemma 24 we have, in polar coordinates, 

z) ~ r^tQtie) + tr^-''hi{t, r, 9) + r^-^Rk{9) + r^-^^(r, 9)), 

with R G C~(E+ X §2n-i) and hi G C°°(M x M+ x S^^-^). We choose new 
coordinates 

Xo(i, r, 9) = igt(^) + r'^-^ithdt. r. + Rk{9) + r^(r, ^)), 

(Xl,X2,-,X2n)(^,?^,^) = (?^,^')- 

Hence I D(t'r e) K*-*' ^- ^0) ^ IQo(6'o)| 7^ and in the new system of coordinates 
we have ^(t, z) = (xoXi)(^> ^1 ^) ^ required. ■ 

Lemma 28 // 9q G Cq^ with 9o ^ Cr^. i(/iere exists a system of local coordi- 
nates X, near {t, r, 9) — (0, 0, ^0), such that ^ ~ XoX2Xi='=Xi ^ neighborhood 
o/(xo,Xi,X2) = (0,0,0). 

Proof. As in Lemma 27 we write 

*(t, ;2) = tr2(g,(^) + r^-2/ii(t, r, ^)) + r'{Rk{9) + rR{r, 9)). 

Since Rk{9o) 7^ in a suitable neighborhood of (6*0,0) we have Rk{9) -\- 
rR{r, 9) 7^ 0. Up to a permutation of the variables 9i, we can suppose that 
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^^(^o) 7^ 0. In a neighborhood of {t, r, 9) — (0, 0, ^o), we choose coordinates 



Xo{t,r,e)^t\Ru{e) + rR{r,e)\-, 
Xiit,r,0)^r\Rkie)+rR{r,e)\i, 
X2{t,r,e) ^ Qt{e) +r''-'hi{t,r,e), 

{X3,-,X2n){t,r,e) = (^2,.-,^2n-l)- 

The corresponding Jacobian, at the point (0, 0, ^o), is 



Dx 



(0,0,^o)H|i?ik(^o) 



(^o)|7^0. 



(39) 



In these new coordinates the phase is ^'(t, r, 0) = (X0X2X1 + Xi){'t, 



Lemma 29 If 6q G Cq^ H Cr^, and if'VQo{9o),'VRk{0o) are linearly indepen- 
dent there exists a system of local coordinates x near {t,r,6) = (0,0. ^0) such 
that ^' ~ XoX?X2 ± X1X3 in a neighborhood of (xo, Xi, X2, Xa) = (0, 0, 0, 0). 

Proof. Near ^0 we can complete {Qo{9), Rk{0)) to a system of coordinates on 
the sphere. Up to a permutation of the 6i, we can choose 

(Xo,Xi)(^,'^,^) = {t,r), 

X2{t,r,9)^Qt{9)+rhi{t,r,9), 

X3it,r,0)^Rkie)+rR{r,0), 

(X4, X2n)(t, r, 9) = (^3, -, 02n-\)- 

Then the corresponding Jacobian is 



'D{t,r,e) 



I (0,0,^0) = 



/dx2 0x2 \ 

dOi 002 
9X3 9x3 



(0,0,^0) = 



a^i dd2 

dRk dRk 



h) ^ 0. 



In these new coordinates we have ^'(t, z) = (X0X1X2 + XiX3)(^) ' 

Remark 30 Coordinates x form a system of admissible charts near {T, zq) 
and these are singular in z = Zq as coordinates on T*W^. In the three systems 
of coordinates the measures are r'^'^~^\j^j^^^(t,r,9)\dtdrd9, this term r^""-^ 
plays a major role since the critical sets of our normal forms are {r = 0}. 

Combining Lemma 27, 28 and 29 gives 

Theorem 31 If (H^) and conditions of Lemma 29 are satisfied and if T is a 
total period of d^ti^o), the phase function S{t,x,^) — (x,^) + tE^ has one of 
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the following normal forms on the blow-up of {T,xo,Co) '■ 

first normal forms : (±xoXi) '^^O'f {U\Cqq), (40) 
second normal forms : {xoxlx2 ± Xi) near (Cq^XCrJ, (41) 
third normal forms : {xoxlx2 ± XiXa) near Cq^ n Cr^. (42) 

We end this section with two lemmas on asymptotics of oscillatory integrals. 
Lemma 32 There is a sequence {cj)j e X>'(]R+ x IR) such that for k > 1 

j ( j e'^'^-'ait, r)dt)dr ~ ^""^^,(0), (43) 



j=0 



where : 

k l\ 

l+l-k 

with Ai = T{x_ * ) and X- = max(— x, 0). 



ci-^-^l{m®5'i\r)). (44) 



Proof. We define ^(r, r) = Tt{a{t, r))(T), where Tt is the partial Fourier trans- 
form with respect to t. Then we obtain 

00 00 00 

j ( j e'^^'\{t, r)dt)dr = j g{-Xr'', r)dr = A'^ j g{-r^ , ^)dr. (45) 

M '^'^ 

Taking the Taylor series in r of ^(r, r), at the origin, gives 



^(-^'' Tt) = E 0) + A-^i?(r, A), 
Straightforward computations shows that 

,_i /■ N , I ^ A ^ k ' f (j'-g . . i+i-k N+i 

" J 'Y^^ k^~ir~ J ' dr + 0{X- ). 



1=0 



l + l-k 

With x_ = max(— a;, 0) and A;(r) = J^(a;_ * )(r) the lemma follows. ■ 

Lemma 33 There is a sequence {cj)j e ©'(M^ x M) such that for k > 1 

°? 00 

/ e'^'''a{r)dr - ^ A-'^^c.^a), (46) 
J=o 

with : 
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Proof. We use the Berstein-Sato polynomial, see e.g. [16]. We write 

oo oo 
7 

with 7 —]c — ioo, c + ioo[, Re(c) < k~^. Since for all positive r we have 

ak k 

^ki^'y^'-'-'^uu-kz), 

we can compute the asymptotic by the residue method. All poles are simple 
and, by pushing of the complex path of integration to the right, we obtain 



n (j - kz) i 



Straightforward computations then show that /x; = (— l)'^r(^) exp(i7r^). 



6 Proofs of the main theorems 



We start with the simpler case of T being a total period of d^ti^z^). Afterwards 
we study the contributions of non-total periods T, distinguishing out particular 
case of a function 'p^ whose restriction to the linear subspace has constant 
sign. In the following we suppose, without loss of generality, that the support 
of the amplitude contains only one non-zero period of the linearized flow. 



6.1 Blow-up and partition of the sphere 



We will apply the results of sections 4 and 5, and we recall that 

R{x,O^S{T,x,0-{x,0, 

{t - T)g{t, X, + tE, = {t- T){Qt{x, + h{t, x, 0). 

By a time translation and with the polar coordinates (x, ^) = r9, 9 E E>' 
we obtain for the top order part of (27r/i)"72(-E'c, h) : 

I{T,h)^ I j a{t,r9)e^''^'''^^\^''-^dtdrde, 

Mx[0,Oo[S2n-l 

where d9 is the standard surface measure on the sphere. 
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Peirtition of unity on the sphere. 

Since Cq^ fl S^"~^(R) is compact, we can introduce a finite partition of unity 
Y^ilW) + E^K^) +E^fW = on §^«-i(M), 

iG/ j&J l€L 

witli the property that Qt{6) ^ onlj supp(rii^), Rk{6) 7^ on Usupp(l^j^) 
and Cq^nCRf, C U supp(r2f). We split up the integral /(T, h) according to this 
partition of unity and use the normal forms of Theorem 31. On any chart let 
be Jx the Jacobian of the relevant diffeomorphism of blow-up from Theorem 
31. For j in /, J and L respectively, we define 

aiAxo,Xi,-,X2n) = {x''rmO)a{t,rey-'\Jxr'), i e {1,2,3}, (48) 

Then, from Theorem 31, we obtain the local contributions 

J n] {0)a{t, re)ei'^^'''-%^''-^dtdrd0 
= I e^^^^'^"A,,i{xo:Xi)dXodxi: 

IRx[0,oo[ 

for the first normal forms. Also 



J n]{e)a{t, re)e^'^^'''^%^''-^dtdrde 



SxR+xS2"-i 



= / ^2j(xo, Xl, X2)e^(>^«?^^±^^)dxo^^XirfX2, 

KxR+xM 

for the second normal forms and for the third normal forms 

J nf{0)a{t, re)ei'^^'''-%^''-^dtdrd0 

Mx]R+xS2"-i 

= / ^3,Kxo, ••,X3)e^(^°^?>^^±^'^^)rfxo..ciX3, 

MxIR+xR2 

where the amplitudes are respectively given by 

^l,i(X0, Xl) = j «2,i(X0, Xl, X2n)dX2---dX2n, (49) 

Aj(X0, Xl, X2) = j a2j(X0, Xl, -, X2n)dX3-dX2n: (50) 
^3,/(Xo, Xi,X2,X3) = a^AXo^Xi, ■■■,X2n)dXA---dX2n- (51) 



It is convenient, for the calculations below, to introduce Aj^i = Xi^^ 

for j e {1,2,3}, c.f. Remark [30]. By construction the functions j are of 

compact support in their system of coordinates. 
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Remcirk 34 We obtain for each new phases the following critical sets 

^iXoxl) = {Xi = 0}, 
' ^(XoX?X2 + X?) = {xi = 0}, 
^(XoX?X2 + XiXs) = {Xi = 0}. 

v 

Where €{f) denotes the critical set of a function /. 



6.2 Analysis in the case of a total period 



First normal forms. 

We note and J-'t the total and partial Fourier transform with respect to t. 
For an amplitude of the form a{r,t) = C(r^"~^) and k — 2, where r — xi, 
Lemma 32 shows that the first non-zero coefficient is obtained for Iq — 2n — l. 
With A = h"^, this gives the contribution 

A~" 

^^^_-^^, C2n-l(A») = (A2n-1 0)) , 

where h.2n-i{t) — T{x^^){t). If we define the distribution 

= (A2n-l(Xo),y'(x"')*(^^l(^)«(^,^^)l>^Xr')(XO,0,X2,...,X2n)c?X2-C?X2n), 

we obtain, as A — > +cxd, the asymptotic equivalent 



oo 

e 







*'^"^'^i,i(Xo, Xi)dXodxi = -^Ai,i(a)A-" + 0(A-"-^). (52) 



Second normal forms. 

If we write a for yl2j, A = and (t,r, f) for (xO)Xi)X2), then we have to 
analyze the asymptotic behavior of the oscillatory integral 

oo 

h{X)=j J e'^^'''''^''''^a{t,r,v)dtdvr^''-'dr (53) 

R2 

If we let a(r, r, f) = jF((a(t, r, f ))(r), the Fourier transform with respect to t, 
then we find by easy manipulations, that 

oo 

/^(A) = A-i J r^^-^e'^'-' J d{-w, r, ■^)dwdr, (54) 
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where we made the change of variables {v, r) — > {w, r), w — Xr^v. 
We now Taylor expand up till order n — 3 : 

w w 
^(~'^' ^) ^ 12 ^ji-'^^ ^-'^ Ar^-*^ ^ -Rn-2(r, w, A), 

where 



and 



Substitution in Eq.(54) leads to 



where 



11—3 
3=0 







and 



oo 

r„_2(A) = A-^ 1 1 Rn-2{r, w, X)e'^'\-^''-''drdw. 

M 

Observe that the oscillatory integrals Jj(A) can be treated individually using 
Lemma 33, but we first analyze the remainder term, r„_2(A), which equals 

OO 1 / \ Q 



A-1 / ''"^ 1 1 -""^""«(--, r, 0) V^iyr^^^-^- (55) 

R ^ ' 

We split the integral with respect to dr as : 

oo yl oo 

dr = j dr + j dr, 

A 

where A = A{X) will be chosen below. We accordingly split r„_2(A) as ?"rj'^2(A) + 
r^'^2(A) and r^'^2(A) is given by Eq.(55), with the integral restrained to [0, A]. 
Easy estimates then show that r^'_2(A) < CA^A"*^""^^ with C independent of 
A. 
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Next, for rn-2{^) we do an integration by part with respect to t, this leads to 

°° n-2 



A 



/ -^11 w-^d:-^ai-w,r, ^Y-^^A^dtdwdr. (56) 

A ]R ^ ' 

We then observe that the first integral in Eq.(56) is equal to 
7 J (n — 2)! 

]R ^ ^ 

= J„_2(A) + C»(A2A-(--i)), 

by similar estimates as for r^'^2(A)- Finally, the last integral in Eq.(56) can be 
estimated by 



A f r \w 



(n-2)! 



n-l 

— ||9r'a(-«;,r,.)||oo«r < CA-"|log(A)|, 



A M 

remembering that a has a compact support in In conclusion, we find that 

n-2 n-2 

^2(A) = E J,(A)+C»(A-("-^)A2)+C»(A--log(^)) = ^ J,(A)+C»(A-"log(A)), 
where we have chosen A — X"^. 

Now by Lemma 32, and since ao(r) = yl2j(0, r, 0), we obtain 

J„(A) = ^.oA^^(0) + 0(A^), (67) 
with jiQ — — |r(^^)exp(i7r^^). While, if j > 0, the same lemma shows that 

2n + (j + l)(fc-2) 

J,(A) = 0{\ ^ ), < J < n - 2. 

All of the latter are dominated by the top order term of Jo(A), and dominate 
the remainder, which is (!?(A~"^log(A)), since for j = n — 2 and for all /c > 3 : 

2n + (/,•- 2) (j + 1) ^ 2 

^ — L = n-l + -<n. 

k k 

As concerns the error term in (56), this also dominates the remainder term, 

since n > for all n > 2. Since = A^^. wc have shown that the contribution 
fe— 1 — 

of the second normal form to the asymptotics of /(T, h) is : 

h{h) = /io^,i(0, 0, 0)h—^ + 0{h^^). (58) 
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Third normal forms. 

We use the same strategy as for second normal forms. If A = h~^, {t, r, v, s) — 
(Xo,Xi,X2,X3) and a = igj, we write 



oo 



/3(A) = /(/ e'^^''^'''+''''^a{t,r,v,s)dtdvdsy-^dr, 

M3 

with a(r, r, v, s) — J^t(a(i, r, v, s))(t), a Taylor expansion, up till order n — 3, 
gives again 

w w 
«(-w,r, — ,s) = Y.^Cj{-w,r,s){—y + Rn-2{r,w, X, s), 

with 

Rn-2{r,wA,s) = ^r-M-y^,r,^^,s)^j^^dt, 

^ ^' 

Cj{w,r,s) = -^(w,r,0,s), 

f ■ d"^ ■ 

=y Cj{-w,r,s)w^dw^ {-iy{{-^^ya){0,r,0,s). 

This leads to 

n— 3 

73(A) = E^i(A)+rn-2(A) 

i=o 

n— 3 „ 

= ^ A-(^+^'^ / e'^'^''aj{r,sy''-^-^^drds + rn-2{X)- 
j=o ■' 

The Kj{\) can be treated individually via Lemma 32. For j — 0, we obtain 
Ko{X) = A-^y^ e^^"'*ao(r, s^-^drds 

M 

~ A-^(f: A-^Q(ao(r, s)r^--=^) + C»(A-'^)). 

The leading term, obtained for I — 2n — 3, is 

A-(^)c2„_3(ao(r,s)r2«-3) = --A-(^)^.F(x_ ao(0, s)^ . 

Hence for our amplitude we have the main contribution 

-Ia-(^) (.F(xr^)(x3), 4^0, 0, 0, X3)) + 0(A-(^)). (59) 
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Like for the second normal forms, the other terms Kj[X), with j > 0, and the 
remainder r„_2(A) give contributions of strictly lower orders. 

Finally, on each local chart the main contributions are 

first normal forms : co,i(a)A~", 
< second normal forms : 00,2(0) A ^ 
third normal forms : 00,3(0) A 

The contributions of charts 2 and 3 are dominant, since > — n, V/c > 2. 

Remark 35 The proofs above show in fact much more than just an asymp- 
totic equivalent for I{T,h), and therefore for 72(£'c, /i). They show the exis- 
tence of a limited asymptotic expansion in the case of indefinite Qt 

I{T,h) =Y.c,h'''^^ +0{h-\log{h)\), (60) 

where the sum is over all u such that '^'^+'^^^+'^ < n, or u < (A; — 2) (n — 1), and 
of a complete asymptotic expansion if Qt is definite. A similar remark applies 
for the case of a non-total period, which we examine in the next section. 

We now compute the leading term of the expansion in case of a non-definite 
Qt and for T a total period of the linearized fiow. 

Case of an empty intersection of cones. 

With Rk on Cp2\{0}, we can assume that Rk is positive on Cp^\{0}. The 
main contribution is here given by the second normal form and is 

/ A,jixo, XI, X2)e^'(>^°^?^^+>^^)dxorfXirfX2 

Rx[0,oo[xR 

, 2-fc-2n ~ . . „ , 2-fc-2n 1 . 

= //feA^^A2,,(0, 0, 0) + 0{X— -k), 

with Ilk given by Lemma 33. By definition of %, the amplitude ^2j(0, 0, 0) is 

|(x-')*(^l?(^)|i?fe(^)+r^(r,^)r^a(i+r,re)|Jxr^)(0,X3,..,X2n)(iX3..rfX2n. 
For 2; e M^, we write this delta-Dirac distribution as an oscillatory integral 

_1_ J e-^(^'(^0'^i'^2){t,r,(9)>^2^^^)|^^^^^) ^^^^^^^^^ re)dzdtdrde. 

If we use 1/2 = Z2\Rk{0) + ri?(r, 6*)!^, integration in (1/2, t) gives 

(27r)'i2j (0,0,0) = j e-'«"i'^3)'(^°'^^)(*'°'^)>Q5(^)|i?fc(^)r^a(t+r, 0)dtdedzidz3. 
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Since Xo{t,0,9) — t\Rk{0)\ »=, with yi — zi\Rk{0)\ by integration in {t,yi) 

(27r)i2j(0,0,0) = a(T,0) J e-'''P^^'^{e)\Rk{e)\-"-^dedzs. 

We define the Liouville measure dLp^ on CpjHS^""^ via dLp^{6) Adp2{0) ~ dO. 
Since charts associated to second normal forms cover the trace of the cone, by 
summation over the partition of unity we obtain 

2n-\-k 2 f 2n 2 1 

I{T,h) = /i^(/ifca(T,0) J \Rk{e)\-—dLp,{e) + 

Cp2nS2n-l 

where 6 are now local coordinates on the surface CpjlHS^""^. The top-order con- 
tribution to the trace formula follows from 72 (-Ec, T, h) = (27r)~'^~^/i~"/(T, h). 



Case of a non-empty intersection of cones. 

Here the main contribution is given by the normal forms 2 and 3. The local 
contribution of any chart associated to second normal forms can be computed 
like in the previous section. The contribution of the third normal forms is given 
by equation (59) with the amplitude A^ ilO, 0, 0, Xs)- Let be 2; = {zi, Z2, Z3), we 
use again an oscillatory representation of the delta-Dirac distribution via 

_1_ I e-(.,(xo,xi,X2)) (^J^{x'^){x3),As,i{xo,..,X3))dxodxidX2dz. 
Since (xo,Xi) = {t,r), integration w.r.t. {zi, Z3,xo,Xi) gives 
2n{^T{x_ " )(X3),^,K0,0,0,X3)) 

// 2ra — 2 — fc \ 

e-*-2X2(o,o,e) u )(^^(o, 0, 9)), nf{e)a{T, 0)j dedz2. 

A classical result, see [8] volume 1 page 167, is 

2n-2-fc 2n — 2 Tl — 1 2n-2 

^{x- " )(X3) = r(— ^) exp(^7r-^)(x3 + ^0)- — . (61) 

By construction (x2, X3)(0) 0; ^) — {P2{0) , Rk{&)) and with the notations of 
Theorem 4 we obtain, using again p2 as a local coordinate 

/^fea(T,0) j {Rk{e)+i0)-'-^nU0)dLp,{9). 

2n — 2 2n — 2 

But {Rk{9)+i0)-— = Rk{9)-' k on a chart associated to the second normal 
form and the total contribution arising from normal forms 2 and 3 is 

I{T,h)^h'^{l^,a{T,0) J {R,{e)+iQ)-"-^dLp,{e) + 0{h^)). (62) 
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Since a(T, 0) = (p{T) exp{iTp^ (zq)) , see formula (66) below, this proves Theo- 
rem 4 for a total period. ■ 



6.3 Case of a non-total period 

Let be T a non total period of the linearized flow. We can assume, up to a 
permutation of coordinates, that dr = = {zi, Z2) € M'^ x M^"~'^ / Z2 = 0}. 
We can apply the Morse lemma with parameter since the phase function at 
time T is only degenerate in z along (cf- Corollary 10). The quadratic part 
S2{t, X, ^) of the function S is given by 

^tid^s, = d^t{Q){d^S2, + o{\\ (x, 1 T) = (x, dA) + o{\\ (x, 1 T). 

Using Theorem 5.3 of [9] we can assume that det[^^] 7^ and the function 
5*2 is well determined locally. Then we have the following facts 

4^(t, z) = Q <^ Z = Zq, Wt, 

Hess^2(^')(r,2;o) = ^^^^'(T, 2:0), is invertible. 

By the Morse lemma, after a change of variable z — > S and calling z again z, 
we have 

^{t, z) = q{z2) + zi, Z2{t, zi)) = q{z2) + ^1), 
and by Corollary 10 again, q = ^qx = P2\dT- In the following, we note 

R{t,zi) = R{zi,Z2{t,zi)), 
g(t,zi) = g{t,zi,Z2{t,zi)), 

^{t, Zi) = ^{t, Zi, Z2it, Zi)) = R{zi, Z2it, Zi)) + {t- T)g{t, Zi, Z2{t, Zi)). 

With these conventions we can write 

I{T,h)= J J e^«^(^^)a(i,Zi,^2)e^('^(*'"i)+^*-'')^(*'"^»rftd-Zid-Z2, 

where a is the new amplitude after change of variables due to the Morse 
lemma. The stationary phase method applied to the 2;2-integral gives 

r ^ 

J e^'^^^'^^dit, zi, Z2)dz2 = E c./i"+"-''M.(t, ^1) + 0{h''+^-'^+^), 

(63) 
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and in particular for the leading term we have 

^o(t,^i) = a(t,^i,0)|^^|-\ 

d{zi,t) 

(27r)"-'^^exp(ifsgn(gT)) 

I det(5T)|^ 
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We now distinguish two different cases : the quadratic form Qt, that is p2 
restricted to ^t, is definite or non-definite. In the first case only the normal 
forms ±XoXi will occur, i.e. the knowledge of d^ti^o) is sufficient. In the second 
case the main contribution involves Rk, and hence the operator d''~^^t{zo). 

The case Qt definite. 

Up to a complex conjugation we can assume that Qt is positive. With polar 
coordinates Zi = {rO) and a|22=o = 'r6',0), the new amplitude is 

^(Xo,Xi) = / {x-yHt,re,Qy--'\Jx{t,TM)dX2...dx26^. (64) 

with Ao(xo,Xi) = Xi^~^A(Xo,Xi)- Lemma 32 shows that 

oo 

/ / Ao(Xo,Xi)e^^°^^o?Xo^^Xi = --/i'^-(^(xl--^)(xo),Ao(Xo,0)) + 0(/i'^-+^). 

M 

With A(xo) = ^(xl^~^)(xo), substituting the definition of Aq, gives 
;A(Xo),io(xo,0)) = 27r / e''' K{xo{t,r,9))a{t,rd,Q)dtdrdedz. 



But by construction xi = f-i and since we have localized the amplitude near 
T 

(A(xo),io(Xo,0)) = -e'5('^--i)r((iT) / {{tQt{,e)-iQ)-''\a{t + T,Q>))dedt. 

(65) 

We can now use a result of [6], also used in [10] and [2]. Since the propagator 
Exp(|iP/i) is a FIO associated to the Lagrangian manifold of the flow, 

A = {(t,T,a;,d/,r/) / {x,i) = ^t{y,v), ^ = -p{x,0}, 

it's principal symbol in the coordinates {t, y, rj) is given by the half-density 

t 

exp(i j p\^,{y,r]))ds)\dtdydr]\^. 



The representation of the propagator with the kernel 

_J_ J ei(^'~''^'^^-(y'^^\a{t,x,7^) + ha,{t,x,7^,h))dy, 
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leads to the half-density 

\dtdydr]\^ 



a{t,x,r))\dtdxdr)\2 ^ a{t,x,r))- 

I ^^K^x,r])\'^ 

For our unique critical point Zq we obtain 

t 

a{t, zo) = I det(S'l^)(t, ^o)|^ exp(i J p\^s(zo))ds) = exp(itp\zo)) . (66) 



If there is no period of d^ti^o) on supp(<^) Theorem 5.6 of [6] gives 

n 27r J \det{ld- d%{zo)\2 

for a certain mo- The denominator has a zero of order in i = T, hence for 
alH 7^ T in a sufficiently small neighborhood of T, we have 

I det(ld - d^t{zo)\2 
Finally, since the contribution is smooth in t, we obtain that 

r(dr) exp(zfsgn(gr)) dr - I . ,^ 

K{T) = — ,.,,|. - 1 exp(z7r— — sign(gT)), 

A(.)^((.-T-,o)-^ ('-y^;-ff';» ^w\. 

\ |det(Id-d$t(zo)|2 / 

and this completes the proof of Theorem 3. 

Remark 36 If there are no rational relations between the eigenvalues of 
and all contributions of the non-zero periods of d^t{zo) are given by The- 
orem 3, as is shown by Proposition 15. This gives a total contribution 

E l{E,.T,h)= E {K{T)KT{^)-rO{h^\ 
Tesupp(.^)\{o} resupp(<^)\{0} 

where the summation is over the non-zero periods of d^ti^o). 
The case Qt indefinite. 

Here we must study carefully the solutions of the implicit function theorem. 
If Zq is zero then, by uniqueness, we have Z2{t, 0) = for all and 

dz2'^{t,Zi,Z2{t,Zi)) = 0, 

det(HesS;,2*(7',0)) 7^ 0. 
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For our unique critical point zi —0 we obtain 

Equations (22) show that dtZ2{T, 0) = and d,^Z2{T, 0) = 0. Hence near (T, 0) 
we have Z2{t, Zi) — 0{\\{t — T, Zi)\\'^). We use, again, a decomposition w.r.t. T 

^(t, zi, Z2{t, zi)) = zi, Z2(T, zi)) + (t- T)g(t, Zi, Z2(t, Zi)). (67) 
Since we have d^^git, Zi, Z2{t, Zi)) ~ 0, it follows that 

liess,,{g{T, zi, Z2{T, ^i))|zi=o = Qt- 

Proposition 23 gives by identification of homogeneous terms of same degree 

^(T, ^1, Z2{T, zr)) = Rkizi, 0) + 0{\\z,\\''+'). (68) 

The stationary phase method in Z2, see formula (63), shows that 



N 

I{T,h) ~ ^ c,/i"+"-'^'' / A,(t,zi,0)e^(-^(*'"i)+(*-^)^(*'"^»didzi, h^O. 

v=0 •' 



Copying the construction for a total period we obtain normal forms for the 
phase ^{t, zi) = R{t, zi) + {t — T)g{t, zi), with the decomposition w.r.t. Cq^ : 



ylo(i,-2i,0)ei*(*'^i) 



^i,i(Xo,Xi)e^''^°^' outside Cq^, 
^2,,(Xo,Xi,X2)e^(>^°>^?>^^±>^?) near Cq^\Cr^, 
^3,KXo,Xi,X2,X3)e^(^°^?^^±>^^>^3) near Cq, n Cr^. 



Now the dimension is It = 2dT and results obtained for total periods show 
again that the contributions of normal forms 2 and 3 are dominating these of 
normal forms 1. Combining this with the leading term of the stationary phase 
method gives 

I{T, h) = r( / a{t, z^ 0)e*(^(*'^i)+(*-^^^(*'^i»dtdzi+O(/i)), 

I det h 



and the contribution of a non-total period is computed by restriction of all 
objects to '^T- This proves Theorems 3 and 4 in their general forms. 
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7 Examples 



Perturbation of harmonic oscillators. Let be 



(69) 



with critical energies and We consider Eq = {{x, ^) / H{x, = 0}, the 
origin is the only singularity of this surface. The system is integrable with 



xi cos(i) + ^1 sm{t) 
X2 cos((4(xi + ei) - l)t) + 6 sin((4(xi + ^1) - 1)^) 
^1 cos(t) — Xi sm(t) 

\ 6 cos((4(a;i + {|) - l)t) - x^ sm{{A{xl + {|) - l)t) J 



On Eo the function 4(a:|+^|) — 1 is bounded and some elementary arithmetical 
considerations show that wc can find a neighborhood V{T) of the origin, in 
Eo, such that all periodic trajectory has a period greater than |T|. Hence (H^) 
is true here. We have 



( cos(i) sm{t) \ 

cos(t) -sin(t) 

-sin(t) cos(i) 

y sin(t) cos(t) J 

and 27r is a total period. Expanding the flow in a Taylor-series gives 

/ ^ 

46(^1 + el) 







\-^X2{xl + Ci) J 

Hence R^lx,^) is negative on Cq\{0}, since 



Siegel and Moser example. To illustrate some properties on periodic tra- 
jectories, Siegel & Moser introduced in [12] the Hamiltonian 

H{x, = lixl + e,) - {xl + el) + x,C^X2 + l{xl- (70) 
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The origin is the only critical point on Eq = {{x,^) / H{x,^) = 0} and 
trajectories from the surface {xi = = 0} are 27r-periodic. Withp = {xl+C,i), 
q = {xl + ^1), we obtain, p" = ipq + (see [12]). Since p is strictly positive, 
p is strictly convex and hence non-periodic. Consequently, on Eq there is no 
non-trivial periodic trajectories. The linearized flow at the critical point is 



cos(i) sm{t) 



cos(2t) sin(2i) 
-sm{t) cos(i) 
\^ -sin(2i) cos(2i) J 

with a resonance of order 3 since, 2wi —W2 — O. For the "period" 27r we have 



The intersection Cq fl Cr^ is obtained by solving the system 

{xi+ei)-2{xi+e2)^0: 



This leads to the surfaces 



X2{xi,Cl) = 



I xia 



6(2^1,6)^-^2 



72X1^1 V 3^1 + 6 ^ 



X2{xu^l) = 



•^1 SI / •^I'il 



V2xieiM+ 6' 
6(..6) = V2gfJ 



,:^f + 6 

By symmetry we just examine gradients on the first surface {Si), we have 



/ 



VR: 



■3|5i 



TT 

4 



-V2a;i6'Y- 
2xi6 



X 



it 



\ 



vg. 



Si 



I 



\ 



2xi 

^ 4 -a [ 

V2xi6V^? + 6' 
26 

-4\/2. 



xi+ei ) 

But, since the minor determinant extracted from VQi^i and Vi?3|5^ 



i^(^i,6) 



V 



-:^i6'^ 



\ x{ + 6^ 
^?6^ 



^?6 



^?6^ / 



^ x\i\ 



-xMxl+e,): 



is non zero for {xi, 6) 7^ 0, VQ, V-R3 are hnearly independent on CqHCr^ flS^ 
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